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ABSTRACT 

Employing  the  autoregressive  (AR)  technique  and  the  principle  of  displaced  phase 
centre  antenna  (DPCA)  we  construct  an  optimum  adaptive  DPCA  processor  for 
moving  target  detection  from  airborne  phased  arrary  radar  data  collected  under  non- 
DPCA  conditions.  It  is  fundamentally  different  from  the  existing  adaptive  DPCA 
which  is  not  optimum.  The  number  of  range  samples  it  needs  to  estimate  its 
parameters  is  only  approximately  twice  the  number  of  antenna  elements,  significantly 
smaller  than  the  number  required  by  the  conventional  space-time  adaptive  processing 
(STAP)  and  other  algorithms.  Computationally  it  only  requires  5-10%  cost  of  STAP. 
The  processor  is  tested  using  both  the  simulated  and  genuine  airborne  phased  array 
radar  data.  With  ample  samples,  its  performance  approaches  optimum.  In  the  case  of 
reduced  samples,  it  considerably  outperforms  STAP  and  others  examined. 
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Phased  Array  Radar  Data  Processing  Using  Adaptive 
Displaced  Phase  Centre  Antenna  Principle 


EXECUTIVE  SUMMARY 

The  Australian  Defence  Force  (ADF)  is  acquiring  an  Airborne  Early  Warning  and 
Control  (AEW&C)  capability  under  Project  AIR  5077  (Project  Wedgetail).  The  AEW&C 
aircraft  is  equipped  with  an  L-band  phased  array  radar  system  for  moving  target 
detection.  An  electronically  scanned  phased  array  antenna  offers  rugged,  reliable  and 
comprehensive  detection  potential.  Flowever,  how  to  process  massive  phased  array 
radar  data  in  given  time  frames  remains  a  key  issue  in  phased  array  radar  data 
processing. 

Although  space-time  adaptive  processing  (STAP)  has  attracted  a  lot  of  attention  in  the 
research  literature,  its  computational  demands  have  limited  its  application  in  real-time 
airborne  radar  systems.  In  this  report  we  investigate  a  new  technique  for  adaptive 
displaced  phase  centre  antenna  (ADPCA)  clutter  cancellation  which  yields  almost  as 
good  a  performance  as  STAP  at  a  fraction  of  the  computational  cost.  In  regions  of 
limited  data  samples  it  is  superior  to  STAP. 

STAP  yields  the  optimum  signal-to-interference-and-noise  ratio  (SINR)  for  detection  of 
moving  targets  embedded  in  Gaussian  distributed  thermal  noise,  clutter  (echoes  from 
the  Earth's  surface)  and  broadband  noise  jamming.  However,  there  are  two  critical 
issues  in  real-time  implementation  of  STAP  in  airborne  radar  systems.  First  it  needs  a 
large  number  of  range  samples  to  estimate  the  covariance  matrix  of  the  undesired 
signals,  which  may  not  be  met  sometimes.  The  second  and  the  more  critical  issue  is  the 
computational  requirement.  The  dominant  computation  in  STAP  is  the  inversion  of  the 
covariance  matrix  which  is  not  only  computationally  expensive,  but  also  a  nonlinear 
transform  making  parallel  processing  difficult  if  not  impossible.  Numerous  algorithms, 
aimed  at  reducing  computational  cost  and  sample  data  requirement  have  been 
proposed  in  the  past. 

Displaced  phase  centre  antenna  (DPCA)  processing  is  a  technique  for  countering  the 
platform  motion  induced  clutter  spectrum  spreading.  The  basic  concept  is  to  make  the 
antenna  appear  stationary  even  though  the  platform  is  moving  forward  by 
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electronically  shifting  the  receive  aperture  backwards  during  the  operation.  Since  its 
origins  in  the  1950s  at  General  Electric  where  it  was  applied  to  airborne  early  warning 
(AEW)  radars,  various  DPCA  techniques  have  been  developed.  For  airborne  radars  a 
typical  way  to  achieve  this  is  to  adjust  the  radar  pulse  repetition  frequency  (PRF) 
according  to  the  platform  velocity  so  that  the  first,  second,  etc.  antenna  elements  at  the 
current  pulse  effectively  move  to,  respectively,  the  exact  positions  of  the  second,  third, 
etc.  antenna  elements  at  the  previous  pulse,  and  so  on.  This  is  often  referred  to  as  the 
DPCA  condition.  In  practice,  even  if  the  DPCA  condition  is  satisfied,  the  clutter 
cancellation  is  still  limited,  due  to  various  disturbances  introduced  by  the  radar  (such 
as  imbalance  among  antenna  channels  and  phase  errors,  etc.),  platform  (instability  of 
velocity  and  crabbing,  etc.),  and  clutter  environment  (clutter  intrinsic  motion).  To 
overcome  these  issues,  the  so-called  ADPCA  concept  has  been  introduced.  However  to 
construct  an  ADPCA  processor  for  airborne  phased  array  radars  is  not  as  easy  as  it 
seems.  The  existing  version  of  ADPCA  was  developed  about  ten  years  ago. 
Unfortunately,  we  prove  that  this  existing  ADPCA  is  not  an  optimum  processor  and 
could  lead  to  significant  SINR  loss. 

Inspired  by  the  DPCA  and  parametric  adaptive  matched  filtering  (PAMF)  approaches, 
this  report  derives  another  ADPCA  algorithm  which  is  fundamentally  different  from 
its  predecessor.  Unlike  the  existing  one  which  is  not  optimum,  the  proposed  ADPCA 
has  its  roots  in  the  autoregressive  (AR)  process  and  is  an  optimum  processor.  It  is 
mathematically  similar  to  the  PAMF  approach,  but  has  some  advantages  over  the 
PAMF  processor.  The  approach  does  not  require  the  DPCA  condition  because  its 
parameters  are  adaptively  estimated,  which  automatically  takes  account  of  non-DPCA 
conditions  as  well  as  various  other  disturbances.  It  significantly  reduces  the  size  of 
sample  data  needed  to  estimate  its  parameters  and  the  computational  cost  compared  to 
its  predecessors  and  STAP. 

The  performance  of  the  algorithm  has  been  assessed  using  two  airborne  radar  datasets, 
one  generated  using  the  high  fidelity  airborne  radar  simulation  software,  RLSTAP,  and 
the  other  collected  by  the  Multi-Channel  Airborne  Radar  Measurements  (MCARM) 
system.  These  two  datasets  were  carefully  chosen  to  cover  various  issues.  First  of  all, 
the  datasets  did  not  satisfy  the  DPCA  condition,  so  that  the  adaptability  of  the 
proposed  ADPCA  could  be  examined.  Secondly  the  radar  did  not  look  in  the  broadside 
direction  in  the  RLSTAP  dataset.  The  second  dataset  was  collected  from  a  flight  trial  of 
MCARM  and  included  the  effects  of  aircraft  crabbing  motion  (the  crab  angle  is  as  large 
as  7  degrees).  Other  decorrelation  effects,  including  radar  instability,  clutter  intrinsic 
motion,  range  foldover  and  interference  caused  by  the  aircraft  etc.,  were  automatically 
included  in  the  dataset.  Therefore  the  evaluation  of  the  proposed  processor  is  realistic. 
The  results  of  STAP  have  served  as  benchmarks  in  evaluation. 

It  has  been  found  that  the  proposed  ADPCA  performs  nearly  as  well  as  STAP, 
suffering  at  most  a  few  dB  of  processing  gain  loss  in  the  vicinity  of  the  Doppler  of  the 
mainlobe  clutter.  However  if  there  are  insufficient  clutter  samples  to  accurately 
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estimate  the  covariance  matrix,  the  performance  of  STAP  is  severely  degraded  whereas 
the  proposed  ADPCA  still  performs  as  well  as  before.  Mathematically  STAP  requires 
an  estimation  of  its  covariance  matrix  whose  size  is  the  product  of  the  number  of 
antenna  elements  and  the  number  of  pulses  in  a  CPI.  On  the  other  hand,  the  proposed 
ADPCA  estimates  its  parameters  by  the  use  of  covariance  matrices  whose  size  is  no 
more  than  the  number  of  antenna  elements.  In  addition,  in  estimating  parameters, 
STAP  uses  averaging  processing  only  in  the  fast-time  domain,  whereas  the  proposed 
ADPCA  utilises  averaging  processing  in  both  the  fast-time  and  the  low-time  domains. 
These  two  differences  make  the  proposed  ADPCA  more  robust  and  require  much  less 
sample  data.  In  general,  the  parameters  of  ADPCA  can  be  satisfactorily  estimated  once 
the  number  of  range  samples  is  equal  to  or  greater  than  twice  of  the  antenna  elements. 
The  existing  ADPCA,  on  the  other  hand,  is  not  an  optimum  processor  and  its 
performance  usually  is  poor  and  suffers  significantly  especially  when  the  target's 
Doppler  is  close  to  that  of  the  mainlobe  clutter. 

The  computational  cost  of  the  proposed  ADPCA  has  been  estimated.  In  general,  it  only 
has  approximately  5-10%  of  the  computational  cost  of  STAP.  Since  most  of  the 
computation  in  the  ADPCA  algorithm  is  linear  transforms,  parallel  processing  and /  or 
hardware  realisation  can  be  easily  implemented.  In  contrast,  the  dominant  calculation 
of  the  STAP  algorithm  is  the  inversion  of  the  covariance  matrix  which  is  not  a  linear 
transform  and  limits  the  application  of  parallel  processing.  In  this  sense,  the 
computational  savings  of  the  ADPCA  algorithm  is  even  greater  than  the  simple 
estimation  of  operational  counts  presented  in  the  report. 

In  conclusion,  the  proposed  ADPCA  algorithm  has  significant  advantages  in  requiring 
far  fewer  samples  and  much  less  computation  while  achieving  robust  performance 
with  the  SINR  improvement  close  to  optimal.  It  therefore  has  a  great  potential  to  be 
implemented  in  real-time  airborne  radar  systems. 
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1.  Introduction 

The  Australian  Defence  Force  (ADF)  is  acquiring  an  Airborne  Early  Warning  and 
Control  (AEW&C)  capability  under  Project  AIR  5077  (Project  Wedgetail).  The  AEW&C 
aircraft  is  equipped  with  an  L-band  phased  array  radar  system  for  moving  target 
detection.  An  electronically  scanned  phased  array  antenna  offers  rugged,  reliable  and 
comprehensive  detection  potential.  However,  how  to  process  massive  phased  array 
radar  data  in  given  time  frames  remains  a  key  issue  in  phased  array  radar  data 
processing. 

Although  space-time  adaptive  processing  (STAP)  has  attracted  a  lot  of  attention  in  the 
research  literature,  its  computational  demands  have  limited  its  application  in  real-time 
airborne  radar  systems.  In  this  report  we  investigate  a  new  technique  for  adaptive 
displaced  phase  centre  antenna  (ADPCA)  clutter  cancellation  which  yields  almost  as 
good  a  performance  as  STAP  at  a  fraction  of  the  computational  cost.  In  regions  of 
limited  data  samples  it  is  superior  to  STAP. 

This  Section  briefly  reviews  techniques  and  problems  in  this  area  in  connection  with 
the  technique  to  be  introduced. 

STAP  yields  the  optimum  signal-to-interference-and-noise  ratio  (SINR)  for  detection  of 
moving  targets  embedded  in  Gaussian  distributed  thermal  noise,  clutter  (echoes  from 
the  Earth's  surface)  and  broadband  noise  jamming  (Ward,  1994,  Klemm,  2002). 
However,  there  are  two  critical  issues  in  the  real-time  implementation  of  STAP  in 
airborne  radar  systems.  First  for  a  range  cell  (or  range  cells)  under  test  (CUT),  data 
from  a  large  number  of  adjacent  range  cells  are  required  in  order  to  compute  the 
covariance  matrix  of  the  undesired  signals.  The  available  sample  cells  are  sometimes 
fewer  than  required  for  accurate  estimation  of  the  covariance  matrix.  This,  however, 
may  be  partially  solved  by  assuming  the  structure  of  the  covariance  matrix  to  be 
known  (Steiner  and  Gerlach,  1998,  2000,  Gerlach  and  Picciolo,  2003,  Bresler,  1988).  The 
second  and  more  critical  issue  is  the  computational  requirement.  The  dominant 
computation  in  STAP  is  the  inversion  of  the  covariance  matrix,  which  is  a  nonlinear 
transform  making  parallel  processing  difficult  if  not  impossible.  Numerous  algorithms, 
aimed  at  reducing  the  requirement  of  either  sample  data  or  computational  cost  or  both 
have  been  proposed  (Klemm,  2004,  Guerci,  2003,  Wang  and  Cai,  1994,  Ward  and 
Kogon,  2004,  Wang  et  al,  2003).  Most  of  these  algorithms,  however,  are  based  on 
various  assumptions,  hence  they  are  confined  in  a  sub-space  and  only  partially- 
adaptive  (Ward,  1994). 

Applying  space-time  non-adaptive  processing  techniques  is  another  alternative  to 
minimise  real-time  computational  demand.  Noticing  the  natural  characteristics  of  the 
clutter  ridge  lying  across  the  Doppler-azimuth  plane.  Farina  et  al  (1998)  have 
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constructed  a  nonlinear,  non-adaptive  space-time  processor  which  incorporates  a  set  of 
filters  with  different  notch  widths  corresponding  to  different  space-time  correlation 
coefficients.  In  application,  input  data  are  passed  through  the  set  of  filters  and  the  filter 
with  the  minimum  estimated  clutter  output  power  is  selected  as  the  best  filter  for  the 
particular  input  data.  Using  the  displaced  phase  centre  antenna  (DPCA)  principle, 
Richardson  (1994)  has  shown  that  the  inverse  of  the  covariance  matrix  (ICM)  has  an 
asymptotic  solution,  and  consequently  constructed  an  optimum  weighting  vector  for 
an  airborne  radar  system  satisfying  the  DPCA  condition.  Dong  (2005)  has  proven  that 
the  ICM  is  approximately  independent  of  the  clutter  environment.  As  a  result,  a  pre¬ 
built  space-time  non-adaptive  processor  has  been  proposed. 

Employing  the  autoregressive  (AR)  technique  (Marple,  Jr,  1987),  Roman  et  al  (2000) 
have  adapted  the  parametric  filtering  concept  to  processing  phased  array  radar  data, 
and  consequently  proposed  an  approach,  called  parametric  adaptive  matched  filtering 
(PAMF)  (Roman  et  al,  2000,  Michels  et  al,  2003,  Rangaswamy  and  Michels,  1997, 
Rangaswamy  et  al,  1995,  Robey  et  al,  1992,  Roman  et  al,  1997).  The  key  concept  in  this 
approach  is  that  measurements  of  pulse  m  are  estimated  using  previous 
measurements  of  pulses  m  - 1,  in  -  2,  •  •  • ,  m  -  p  ( p  is  also  referred  to  as  the  order  of  the 
filter,  p  >  1).  The  residue,  i.e.,  the  difference  between  the  actual  measurement  of  pulse 
m  and  the  estimated  measurement  from  previous  measurements  of  pulses 
2, is  then  processed.  Since  the  residue  is  uncorrelated,  its 
covariance  matrix  becomes  a  block  diagonal  matrix,  which  in  turn  significantly 
simplifies  the  processing  and  reduces  the  requirement  of  computation.  This  technique 
has  been  reviewed  in  detail  and  further  developed  by  Dong  (2006). 

DPCA  processing  is  a  technique  for  countering  the  platform  motion  induced  clutter 
spectrum  spreading.  The  basic  concept  is  to  make  the  antenna  appear  stationary  even 
though  the  platform  is  moving  forward  by  electronically  shifting  the  receive  aperture 
backwards  during  the  operation  (Morris  and  Harkness,  1996).  Since  DPCA  had  its 
origins  in  the  1950s  at  General  Electric  and  was  applied  to  airborne  early  warning 
(AEW)  radars,  various  DPCA  techniques  have  been  developed  (Muehe  and  Labitt, 
2000,  Nohara,  1995).  A  typical  way  to  achieve  this  is  to  adjust  the  radar  pulse  repetition 
frequency  (PRF)  according  to  the  platform  velocity  so  that  the  first,  second  etc.  antenna 
elements  at  pulse  m  effectively  move  to  the  respective  positions  of  the  second,  third 
etc.  antenna  elements  at  pulse  m  —  1 ,  and  so  on.  This  is  often  referred  to  as  the  DPCA 
condition.  In  practice,  even  if  the  DPCA  condition  is  satisfied,  the  clutter  cancellation  is 
still  limited,  due  to  various  disturbances  introduced  by  the  radar  (such  as  imbalance 
among  antenna  channels  and  phase  errors,  etc.),  platform  (instability  of  velocity  and 
crabbing,  etc.),  and  clutter  environment  (clutter  intrinsic  motion).  To  overcome  these 
issues,  the  so-called  ADPCA  algorithm  has  been  introduced  in  airborne  side-looking 
linear  array  radars  (Blum  et  al,  1996,  Klemm  2002,  Chapter  7,  Guerci,  2003,  Chapter  5). 
Unfortunately,  we  show  that  it  is  not  an  optimum  processor  in  Section  3  and 
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consequently  could  have  significant  SINR  loss  in  Section  5.  We  refer  to  this  ADPCA  as 
Blum's  ADPCA  in  this  report. 

Inspired  by  the  DPCA  and  PAMF  approaches,  this  report  proposes  another  version  of 
ADPCA  algorithm  applicable  to  airborne  side-looking  phased  array  radar  data 
processing  for  moving  target  detection.  Although  the  name  of  ADPCA  has  frequently 
appeared  in  the  literature,  the  technique  developed  here  is  fundamentally  different 
from  the  existing  ones.  Unlike  the  Blum  ADPCA  which  is  not  optimum,  the  ADPCA  to 
be  introduced  here  has  its  roots  in  the  autoregressive  (AR)  process  (Marple,  Jr,  1987), 
and  is  mathematically  similar  to  the  PAMF  approaches,  but  has  some  advantages  over 
the  PAMF  processor  (refer  to  Section  5  for  the  comparison).  The  approach  does  not 
require  the  DPCA  condition,  because  its  parameters  are  adaptively  estimated,  which 
automatically  takes  account  of  non-DPCA  conditions  as  well  as  various  other 
disturbances  induced.  It  significantly  reduces  the  size  of  sample  data  needed  to 
estimate  its  parameters  and  the  computational  cost  compared  to  its  predecessors.  There 
are  no  assumptions  about  the  form  of  the  clutter  in  the  algorithm,  so  it  is  fully- 
adaptive.  Most  importantly  the  SINR  of  the  processor  approaches  the  optimum. 

Section  2  briefly  reviews  STAP.  When  sufficient  samples  are  available  to  accurately 
estimate  the  covariance  matrix,  STAP  yields  the  optimum  SINR,  and  thus  serves  as  a 
benchmark  for  evaluating  the  results  of  the  proposed  algorithm.  Following  the 
introduction  of  DPCA  and  Blum's  ADPCA  in  Section  3,  Section  4  formulates  the 
proposed  ADPCA  algorithm.  Its  advantages  including  the  requirement  of  much  fewer 
samples  in  estimating  its  parameters  as  well  as  the  computational  savings  are  also 
explained  and  analysed  in  Section  4.  Evaluation  of  the  algorithm  is  presented  in 
Section  5.  Two  datasets,  one  generated  using  the  high  fidelity  airborne  radar 
simulation  software,  RLSTAP,  and  the  other  collected  by  the  Multi-Channel  Airborne 
Radar  Measurements  (MCARM)  system  are  used  in  the  assessment,  with  the  results  of 
STAP  as  benchmarks.  The  results  of  the  Blum  ADPCA  and  traditional  DPCA  with  non¬ 
adaptability  are  also  shown  for  comparison.  Finally  Section  6  concludes  the  report. 

2.  Space-Time  Adaptive  Processing 

This  section  briefly  reviews  the  STAP  principle  as  this  helps  understanding  the 
proposed  ADPCA  later. 

The  mathematical  notation  used  in  this  report  follows  the  convention:  boldface 
uppercase  and  lowercase  letters  represent  matrices  and  vectors,  respectively, 
superscripts  T ,  H  and  *  denote  transpose,  Flermitian  transpose  and  complex 
conjugate,  respectively,  and  C  and  !R  stand  for  the  complex  and  real  number  fields, 
respectively.  For  instance,  x  e  <CA/A  1  denotes  x  to  be  an  MW -element  complex-valued 
column  vector,  and  A  e  (CMxMV  denotes  A  to  be  an  MN  x  MN  -element  complex¬ 
valued  matrix.  The  symbol  0  denotes  the  Kronecker  matrix  product,  and  E{-} 
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expectation.  The  ensemble  average  symbol  is  an  over-hat  A,  an  indication  of  an 
estimate  of  E{-}  in  calculation  using  finite  samples.  Uppercases  M  and  N  are 
reserved  for  the  number  of  pulses  in  a  coherent  processing  interval  (CPI)  and  the 
number  of  antenna  elements  in  azimuth. 

In  signal  processing,  one  often  has  to  estimate  unknown  parameters.  For  instance,  the 
true  covariance  matrix  R„  which  is  in  practice  unknown,  is  often  estimated  from 
range  samples  by  using  the  maximum  likelihood  method.  The  estimated  covariance 
matrix  is  Ri( .  In  consequence  the  weighting  vector  of  a  processor  is  w  corresponding 
to  Ri( ,  w  corresponding  to  R„ ,  and  so  on.  Mathematically  if  the  true  covariance 
matrix  Ri(  is  known,  the  optimum  weighting  vector  of  STAP  is  wsrAp  =  yR^'e 

whereas  the  estimate  is  w  STAP  = ;/  R  „ 1  e  .  In  order  to  reduce  the  number  of  equations 

and  expressions,  this  report  does  not  generally  discriminate  between  w  and  w,  if 
meanings  in  the  context  are  clear  and  often  the  over-hat  (indicating  an  estimate  rather 
than  the  true  value)  will  be  omitted.  In  general,  when  we  deal  with  mathematical 
derivations,  we  assume  the  true  values  of  parameters  are  known  so  symbols  are 
without  the  over-hat.  On  the  other  hand,  when  we  come  to  calculation  or  realisation, 
any  parameters  estimated  from  using  finite  sample  data  should  have  the  over-hat  to 
indicate  the  approximation  or  realisation.  However  if  there  is  no  confusion,  we  will  not 
repeat  expressions  with  the  over-hat. 

Let  the  snapshot  x  =  [xr(0)  xr(l)  •  •  •  xT  (M  -  l)]r  with  x(m)  e  CAxl  be  the  complex 
valued  sequence  of  the  output  of  a  linear  equispaced  N-element  array  after 
demodulation  and  pulse  compression,  corresponding  to  the  return  from  a  single  range 
cell  over  an  M-pulse  CPI.  If  undesired  signals  xu  are  Gaussian  distributed,  it  is  well 
known  that  the  optimal  (i.e.,  STAP)  weighting  vector  is  (Compton,  Jr,  1988,  Ward, 
1994), 

%=^1e  (1) 

where  y  is  an  arbitrary  scalar,  which,  without  loss  of  generality,  we  will  let  equal  1 
and  hence  omit  in  the  following  equations.  Ru  =  E{xux^}  e  CMjVxMV  is  the  covariance 
matrix  of  the  undesired  signals.  The  undesired  signals  xu  may  consist  of  Gaussian 
distributed  clutter,  thermal  noise  and  possibly  point-source  wideband  noise  jamming. 
The  covariance  matrix  is  generally  unknown  and  often  has  to  be  estimated  by  the  use 
of  range  samples  via  the  so-called  diagonally  loaded  sample  matrix  inversion  (DL-SMI) 
method  (Carlson,  1988,  Ward  and  Kogon,  2004),  as. 
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1  K 

Ru  =  — Zx 
Kk ~i 


H 


kXk 


+  8 1 


(2) 


where  the  subscript  k  denotes  range  bins,  8  is  a  small  value  of  the  order  of  the  system 

thermal  noise  level  and  I  is  the  identity  matrix.  The  vector  eeC MNA  is  the  desired 
signal  spatial-temporal  steering  vector. 

e  =  e(/,)®e(/(s)  (3) 


where  ftd  and  fts  are  the  normalised  target  Doppler  and  spatial  (azimuth) 
frequencies,  respectively.  Vectors  e(/,rf)  and  e(/’v)  are  the  desired  target  temporal 
and  spatial  steering  vectors,  respectively,  defined  as, 

e(/J  =  ^=['  exp(j2^)  •••  <zxp{j2n{M  -  \)ftd  )]7  (4) 

and 

e(/J  =  ^=&  exp(j2^te)  •••  exp(j2x(N -\)fts)]r  (5) 


Using  block  component  notation,  (3)  may  be  written  as. 


e(m)  = 


-j=  exp(/ 2^-  m  ftd  )e(/te ) 


m  =  0,1,---,M  -1 


(6) 


The  output  of  the  optimum  processor  is  the  product  of  the  Hermitian  transpose  of  the 
optimal  weighting  vector  times  the  data  snapshot, 

y  =  w>  (7) 


The  coherent  processing  gain  of  the  STAP  processor  approaches  the  upper  limit, 
101og10(MA)  dB,  for  deterministic  signals  whose  bearings  differ  from  the  jamming 
bearings,  and  whose  Doppler  frequencies  differ  from  the  Doppler  frequency  of  clutter 
in  the  radar  look  direction. 


Given  a  target  signal  ae ,  where  a  is  the  amplitude  of  the  target  signal,  the  SINR  is 
defined  as  the  ratio  of  the  output  target  power  to  the  output  interference  and  noise 
power. 
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SINRopt 


\yt  1 
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| 

(8) 


where  =E{\a  |2}  . 


For  an  unknown  signal  amplitude,  a  detection  test  statistic  (DTS)  of  STAP  has  been 
proposed  as  the  ratio  of  the  output  power  to  the  power  of  interference  and  noise 
(Robey,  et  al,  1992), 


A 


STAP 


Iw"  xl2 

1  ^  opt A  1 

1  w "  X  I2 

wH  R  w 

TT  opt  u  TT  opt 

|w",e| 

(9) 


Comparing  (8)  and  (9)  we  can  seen  that  the  STAP  DTS  is  actually  the  SINR  of  the  signal 
x .  In  operations,  once  Astap  exceeds  a  threshold,  target  presence  is  declared. 

3.  Classical  Displaced  Phase  Centre  Antenna 

3.1  Displaced  Phase  Centre  Antenna 

DPCA  processing  is  a  technique  for  countering  the  platform  motion  induced  clutter 
spectrum  spreading.  The  basic  concept  is  to  make  the  antenna  appear  stationary  even 
though  the  platform  is  moving  forward  by  electronically  shifting  the  receive  aperture 
backwards  during  the  operation  (Morris  and  Harkness,  1996). 

An  airborne  linear  equispaced  phased  array  antenna  aligned  with  the  flight  direction 
and  looking  in  the  broadside  direction  is  the  default  configuration  in  this  report 
(though  the  radar  is  allowed  to  steer  its  look  direction  off  the  broadside  direction,  see 
numerical  examples  in  Section  5).  In  particular,  let  a  platform  be  moving  in  the  x- 
direction,  at  a  speed  of  va ,  a  linear  antenna  array  with  N  elements  be  parallel  with  the 
direction  of  motion,  and  looking  at  the  broadside  y  -direction,  as  shown  in  Figure  1. 
Symbols  H ,  R ,  6  and  (f)  denote  platform  height,  range,  elevation  angle  and  azimuth 
angle,  respectively. 
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Figure  1:  Geometry  of  a  linear  airborne  antenna  array. 

With  the  common  full-array  transmit  aperture,  the  clutter  return  from  clutter  patch  P 
received  by  antenna  element  n  for  pulse  in  is  (Ward,  1994), 

cm.n  =  (v2%p)U2  expC/^p ) exp[ /2;r(n  +  2vaTrm / d)^- cos 9 cos <f>]  (10) 

A 

where  cr 2  is  the  thermal  noise  level  of  antenna,  g,>  is  the  single-pulse  single-element 
clutter-to-noise  ratio  (CNR)  for  the  clutter  patch  P ,  cpp  is  a  phase  term  induced  by 
clutter  patch  P ,  Tr  is  the  pulse  repetition  interval  (PRI),  d  and  A,  are  the  antenna 
element  space  and  radar  frequency,  respectively. 

The  DPCA  condition  for  the  above  configuration  requires  the  so-called  clutter  slope 
j3  =  2vaTr/d  (11) 

to  be  an  integer.  When  /?  =  1 ,  we  have  from  (10), 

Cm+k.n-k=Cm,n  k  =  l,2,---  (12) 

It  means  that  the  phase  centre  of  antenna  element  n  —  k  will  move  to  the  phase  centre 
of  antenna  element  n  after  k  pulses  as  shown  in  Figure  2  (a).  As  a  result  the  clutter 
signal  received  by  element  n  from  pulse  m  should  be  statistically  equal  to  the  clutter 
signal  received  by  element  n  —  1  from  pulse  in  + 1 .  The  clutter  cancellation  can  be 
achieved  if  we  subtract  them  from  each  other.  Having  said  so,  it  is  worth  noting  that 
the  nature  of  radar  signals  is  a  random  process.  The  phase  term  cpp  and  the  intensity 
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term  gp  are  random  variables  and  each  obeys  its  own  pdf.  Equation  (12)  therefore 
should  be  understood  as  E{cm+kn_k}  =  E{cmn} ,  or  cm+kn_k-cmn=Smn(k)  where 
8m  n  (k)  is  a  zero  mean  random  variable. 


Rjlse  0 
Rjlse  1 
Rjlse  2 


0 

1 

N-1 

Figure  2:  Effective  array  positions  for  successive  pulses  of  a  CPI  with  /?  =  1 . 


Equation  (11)  imposes  serious  limitations  on  any  implementation  of  the  DPCA 
technique.  Once  the  element  spacing  and  the  pulse  repetition  frequency  (PRF)  are 
fixed,  the  platform  velocity  is  also  constrained.  In  practice,  even  if  the  condition  of  (11) 
is  satisfied,  the  clutter  cancellation  is  still  limited  due  to  other  factors  such  as  irregular 
motion  and  crabbing  of  the  aircraft,  spatial  decorrelation  among  antenna  elements  and 
temporal  decorrelation  among  pulses,  as  well  as  due  to  the  intrinsic  motion  of  clutter. 


3.2  Blum's  Adaptive  DPCA 

Adaptive  processing  can  tolerate  certain  perturbations.  It  is  thus  natural  to  include 
adaptivity  into  DPCA  processing  to  restore  performance  under  non  DPCA  conditions. 
Blum  et  al  (1996)  first  proposed  a  suboptimum  technique  (which  is  optimum  if  the 
number  of  pulses  in  a  CPI  is  two)  and  called  it  adaptive  DPCA.  (Blum  et  al,  1996, 
Klemm,  2002,  Chapter  7,  Guerci,  2003,  Chapter  5).  The  algorithm  was  further  discussed 
and  expressed  in  a  more  general  form  by  Guerci  (2003,  Chapter  5). 

To  derive  the  Blum  ADPCA,  let  us  start  from  the  case  of  a  two-pulse  system. 
According  to  STAP  (1),  the  optimum  weighting  vector  is. 


w 


opt 


=  R0:l 


e(0) 

e(l) 


(13) 


where  R(H  e  (C~,Vx~ v  is  the  covariance  matrix  of  the  undesired  signals  corresponding 

to  pulses  0  and  1.  Vectors  e(0)  and  e(l)eC  ' ' 1  are  sub-spatial-temporal  steering 
vectors  also  corresponding  to  pulses  0  and  1,  respectively,  whose  explicit  expression  is 
given  in  (6).  This  two-pulse  processor  formed  from  the  principle  of  STAP  is  called  two- 
pulse  ADPCA  by  Blum  et  al  (1996).  It  can  be  seen  that  the  Blum  APDCA  processor  is 
the  same  as  the  STAP  processor  and  thus  optimum.  If  the  number  of  pulses  is  greater 
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than  two,  each  pair  of  pulses  forms  a  two-pulse  canceller;  overall  the  weighting  vector 
of  the  Blum  ADPCA  is  (Blum  et  al,  1996,  Guerci,  2003), 


Ro:l 

0 

0 

0 

e(0 : 1) 

0 

Rl:2 

0 

0 

e(l:2) 

0 

0 

0 

0 

0 

0 

D  -1 

1) 

e(M  -  2  ;  M  - 1) 

where  R , 


is  the  covariance  matrix  of  the  undesired  signals  and 


e(m  :  m  + 1)  =  the  spatial-temporal  steering  vector  corresponding  to  pulses 

e(m  + 1)  J 

m  and  m  +  \ ,  for  m  =  0, 2,  •  •  • ,  M  -  2  . 

Equation  (14)  can  be  expressed  as. 


R0:l  0  0 

0  r12  0 

0  0 

0  0  0  R 


where 


0^2Nx2N 

®2  NxN 

®2NxN 

^2  NxN 

l^2Nx2N 

®2NxN 

^2  NxN 

®2  NxN 

M-2^2Nx2N 

where  1 2 Ny2 ,v  a  2 N x  2 N  identity  matrix,  02,Vx;V  a  2 NxN  zero  matrix,  and 


|_e(M  - 1)  J 

is  the  original  steering  vector. 

The  corresponding  test  statistic  is. 
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H 

2 

WADPCAXr 

W  ADPCA^T W  ADPCA 


(18) 


where  xr=Tx  and  R  r  =  TR  ,  r;/ .  If  A  exceeds  a  threshold,  target  presence  is 

declared.  It  should  be  pointed  out  that  because  w  4DPCA  ^  R r'  T e  ( R,1  does  not  exist, 
see  below),  the  denominator  of  (18)  cannot  be  further  simplified  as  in  (9). 

The  above  Blum  ADPCA  processor  (15)  can  be  considered  as  a  linear  transform,  T ,  of 
the  original  vector  space,  i.e.,  the  signal  snapshot  x  becomes  Tx  and  the  spatial- 
temporal  steering  vector  e  becomes  T  e .  However, 


1*0:1  0 

0  r12 

0  0 

0  0 


0 
0 

0  R 


0 
0 
0 

(M-2):(M— 1) 


* 


e{(Tx)(Tx)h  }=  £{xrxf  }=  TRUTH  if  M  >  2  (19) 


Therefore,  the  Blum  ADPCA  is  not  optimal  except  for  M  -  2  . 

The  non-equal  sign  of  (19)  is  easy  to  prove.  One  method  is  given  here.  The  rank  of  the 
matrix  on  the  left  side  of  the  non-equal  sign  is  2N(M  - 1)  as  each  sub-matrix  is  a  full 
rank  2N  x  2N  matrix.  Because  Ru  is  a  positive-definite  Hermitian  matrix,  the  rank  of 
the  matrix  product  on  the  right  side  of  the  non-equal  sign  is 
r(rR»r//)=  r((rR|/2)(rR)(/2)//)=  r(rR^/2)  =  minjr(r),  r(R\!2)}  =  MN,  where  rQ 
is  the  rank  function.  This  proof  also  shows  that  there  does  not  exist  a  general  linear 
transform  T  which  lets  the  equal  sign  hold  except  for  the  case  of  M  =  2  .  The  proof 
also  indicates  the  flaw  in  some  equations  given  by  Guerci  (2003,  Chapter  5). 

To  illustrate  the  above,  we  demonstrate  here  an  explicit  example  oi  N  =  2  and  M  =  3 . 
The  snapshot  vector  and  the  transform  matrix,  respectively,  are. 
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*00 

*01 

*10 

*11 

x20 

*21 


and 


1  0  0  0  0  0 
0  1  0  0  0  0 
0  0  1  0  0  0 
0  0  0  1  0  0 
0  0  1  0  0  0 
0  0  0  1  0  0 
0  0  0  0  1  0 
0  0  0  0  0  1 


(20) 


where  xmn  is  the  echo  received  by  antenna  element  n  for  pulse  in  .  The  transformed 
covariance  matrix  is. 


4r(xx//)r//}=^i 
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(21) 

Since  the  covariance  matrix  is  Hermitian,  we  omit  the  elements  below  the  main 
diagonal  in  (21).  Comparing  (21)  with  (15),  we  can  see  that  the  Blum  ADPCA  assumes 
zero  sub-matrices  (the  unshaded  sub-matrices)  for  all  the  off-main  diagonal  block 
matrices.  Therefore  the  Blum  ADPCA  is  not  optimal  for  M  >  2  .  As  a  consequence,  the 
processing  could  lead  to  a  significant  SINR  loss  if  the  target's  Doppler  is  close  to  that  of 
the  mainlobe  clutter.  We  leave  the  performance  evaluation  to  Section  5. 
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4.  Proposed  Adaptive  Displaced  Phase  Centre 

Antenna 


4.1  Formulation  of  Proposed  ADPCA 

The  starting  point  in  the  proposed  ADPCA  is  to  form  the  difference  between  the 
measurements  of  array  elements  1  to  N  —  1  for  pulse  m  and  array  elements  0  to 
N-2  for  pulse  m  + 1 ,  as 


A  x(m)  = 


Xm,X 

Xm+ 1,0 

A*™,! 

= 

Xm,  2 

- 

Xm+ 1,1 

_^Xm,N-  2_ 

_Xm,N- 1  _ 

Xm+\,N-1  _ 

(22) 


which  can  be  written  in  a  compact  form  as, 

A\(m)  =  [0  I  ¥_|  ]\(m)  -  [i  v  |  0]x(m  +  l)  m  =  0,1,---,M  -2  (23) 

where  OdR1'^1’’1  is  a  zero  vector,  I  A?  l  e  [R( '  1/1  '  1 1  is  an  identity  matrix  and 
Ax(m)e€(A,_1)xl. 

If  the  condition  J3  =  1  is  satisfied  and  there  are  no  decorrelation  effects,  statistically  all 
the  clutter  signals  will  be  suppressed  and  A \(m)  given  in  (23)  will  only  contain  noise 
residue.  If  /?  ^  1 ,  which  is  the  condition  we  are  interest  in,  and /  or  there  exist 
decorrelation  effects,  the  clutter  signal  will  only  be  partially  cancelled.  We  modify  (23) 
as, 

Ax(ra)  =  [0  IA,_1]x(/«)-[A(m)  0]x(m  +  l)  m  =  0,1,--,  A/  -2  (24) 

where  A(m)  e  C|  A  1)x<iV  ^  is  a  parameter  matrix  to  be  determined.  It  means  that 
instead  of  using  a  single  measurement  xm+1  n  ,  to  estimate  xm  n  as  in  (23),  we  now 
using  a  linear  combination  of  xm+l  j ,  •  •  • ,  xm+l  N  2 .  This  is  sometimes  referred  to  as  an 

autoregressive  (AR)  process  (Marple,  Jr,  1987)  or  PAMF  (Roman  et  al,  2000,  Michels  et 
al,  2003,  Rangaswamy  and  Michels,  1997). 

Our  goal  is  to  maximally  suppress  the  clutter,  i.e.,  to  minimise  the  mean  residue  power 
with  respect  to  the  parameter  matrix  A (m) , 
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min  e\axh (m)Ax(m)\  (25) 

A  (m) 


The  optimal  parameter  matrix  A (m)  may  be  found  from  the  least  squares  method  (the 
Lagrange  method), 

o(zr  ]ax  "  fm)Ax(m)  [)  _  () 
dA(m) 


Since  A(/n) in  (25)  is  independent  of  the  expectation  operation,  the  expectation 
operation  and  the  differentiation  are  interchangeable,  and  hence  (26)  can  be  expressed 
as. 


£fa(Ax//(m)Ax(m))|  =  Q 
i  dA(m)  j 


(27) 


The  derivative  of  a  scalar  function  f(Q)  e  (Clxl  with  respect  to  Q  e  <Cmx"  is  an  mxn 
matrix  (Van  Trees,  2002,  p.  1399), 


df/dqu  •••  df/dqln 
of  /  8q  m\  Oif  !dqmn 


(28) 


where  qik ,  /  =1, — ,  m ,  k  =  1,  •  •  • ,  n ,  is  the  element  of  Q  .  After  tedious  but 

d(Ax//(/n)Ax(m))|  _  . 

dA(m)  j 


straightforward  algebra,  the  equation  satisfying  El 


A”  (m)  = 


E{\xm+ifi\  } 

^  {X  m+l,0X  m+l,N—2  } 

-1 

E{Xm+l,0X*n,l)  ■' 

E{Xm+l,0Xm,N-l} 

E{Xm+l,N-2Xm+l,o} 

E{\Xm+l,N-2  } 

E{xm+l,N-2Xm,l} 

^  {Xm+l,N-2Xm,N-\  } 

=  [Rm+1>  m+l  (0:N-2,0:N-2)Yl  Rm+Ul  (0 :  N  -  2, 1 :  N  - 1) 

m  =  0,1,---,M  -2  (29) 


where  R m  „(ml  \m2,n{  :«2)eCn  IMA  11  is  a  block  submatrix  of  the  covariance  matrix 
R1( .  Subscripts  in  and  n  indicate  the  row  and  column  of  the  block  submatrix  of  R  , 
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and  the  ranges  in  the  brackets  indicate  the  rows  and  columns  of  Rm  n .  The  block 
matrix  form  of  R1(  is  given  by. 


R 


®-0,0 


M- 1 


where  R  „  e  C  A  xN ,  m,  n  =  0,  •  •  • ,  M  - 1 . 


(30) 


For  simplicity,  we  rewrite  (29)  as, 

A"  (m)  =  [R'm+hm+l  ]-*  R'm+l  m  m  =  0,-  •  -,M  -  2  (31) 

Since  the  true  covariance  matrix  is  unknown,  A (m),  m  =  0,---,M  -2,  can  only  be 
estimated  using  sample  data  in  practice.  Assuming  the  radar  system  is  well  calibrated 
and  stable,  and  the  pulse  train  is  uniform,  then  for  fast  calculation,  we  can  assume 
A (m)  to  be  independent  of  pulse  number  m  in  a  CPI.  This  assumption  allows  the 
averaging  processing  to  proceed  using  the  slow-time  (pulse)  samples  after  the 
averaging  processing  using  fast-time  (range)  samples,  and  has  advantages  and 
disadvantages.  If  there  are  sufficient  data  samples  and  each  is  an  identically  and 
independently  distributed  (iid)  sample  of  undesired  signal  identical  to  that  of  CUT, 
then  the  averaging  processing  may  lead  to  some  processing  gain  loss  as  temporal 
correlation  effects  as  well  as  differences  are  being  averaged.  On  the  other  hand,  if  the 
data  samples  are  limited,  the  averaging  process  should  make  the  estimate  of  A (m) 
more  accurate  and  result  in  a  robust  processor.  Therefore,  the  assumption  that  A (m)  is 
independent  of  m  is  not  only  for  the  reason  of  reducing  calculation,  but  more 
importantly  requires  less  sample  data.  The  rationale  of  this  assumption  is  given  below. 
If  clutter  is  modelled  by  the  first-order  clutter  model,  i.e.,  the  model  does  not  take 
account  of  various  spatial  and  temporal  decorrelation  effects,  the  covariance  matrix  has 
a  so-called  Toeplitz-block-Toeplitz  pattern  (Ward,  1994),  then  we  can  immediately 
deduce  all  A (m),  m  =  0,  •  •  • ,  M  -  2  to  be  identical.  Temporal  decorrelation  effects  such 
as  range  walk  and  clutter  intrinsic  motion  are  often  modelled  by  the  covariance  matrix 
tapering  (CMT)  method  which  is  often  assumed  stationary,  i.e.,  the  tapering  coefficient 
is  only  a  function  of  the  interval  between  pulses  (Guerci,  2003,  Chapter  4,  Clemm,  2002, 
Chapter  2).  Under  the  CMT  model  assumption,  A (m)  is  still  independent  of  the 
pulses.  Note  also  that  this  assumption  does  not  require  any  specific  spatial  correlation 
form  as  the  averaging  process  is  only  in  the  temporal  domain,  though  most  models  of 
spatial  correlation  are  also  assumed  stationary  (Guerci,  2003,  Chapter  4). 
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An  averaging  process  in  the  slow-time  domain  results  in  the  parameter  matrix  A(m) 
to  be  independent  of  pulses,  as. 


Ah  = 


M- 2 


M  1  m=0 


Yr' 

rvm+l,m+ 1 


i  M-2 

_LyR' 

Ti  /f  1  /Jl+1  J 

M-1S 


M-2 

m=0 


m+\,m+l 


-ir 


M-2 

vr;, 

m=0 


+1 ,/ 


(32) 


The  averaging  process  of  (32)  means  that  temporal  correlations  such  as  the  effects  of 
range  walk  and  clutter  intrinsic  movement  will  be  averaged  (similar  to  the  CMT 
models).  We  will  give  numerical  comparisons  in  the  next  section  to  justify  such  a  slow¬ 
time  averaging  process. 

It  is  understood  that  if  A  =  I ,  then  the  processor  will  reduce  to  the  traditional  DPCA 
with  non-adaptability. 


The  transform  from  x  to  Ax  can  be  expressed  in  a  compact  form,  as 


Ax  =  Bx 


where  B  e  C(M  l)l  A  1)xMN  and  has  the  form. 


0  I  -A  0  0  0  0 

0  0  I  -A  0  0  0 

0  0  0  I  -A  0 


We  perform  the  same  transform  for  the  steering  vector  e ,  as, 

u  =  Be 


(33) 


(34) 


(35) 


Analogous  to  STAP  processing,  for  snapshot  Ax  and  steering  vector  u,  the  optimal 
weighting  vector  is, 

w  =  yRAu  (36) 


where  RA  =  £|AxAx//  J  is  the  covariance  matrix  of  the  residue  signal  Ax  and  can  be 
estimated  using  range  samples,  by 

R4=4£AitAx?  =BR„B"  (37) 

K  k= l 

Since  Ax  is  a  residue  noise  vector,  its  elements  should  be  mutually  uncorrelated  (the 
process  of  (33)  is  also  referred  to  as  a  whitening  process).  The  covariance  matrix  RA 
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therefore  should  only  be  a  diagonal  matrix,  and  its  inverse  is  simply  the  inverse  of  each 
diagonal  element.  In  reality,  the  non-diagonal  elements  of  RAmay  not  be  zero,  but 
their  values  will  be  smaller  compared  to  the  diagonal  elements,  and  the  inverse  can  be 
found  approximately  by  letting  all  non-diagonal  elements  equal  zero,  i.e., 

Ra  ~[diag(RA)]_1  (38) 

The  corresponding  weighting  vector  (36)  becomes, 

w  =yRAu  (39) 


The  weighting  vector  of  (39)  may  not  be  optimal  but  should  be  close  to  the  optimum. 
The  resultant  test  statistic  is. 


A 


ADPCA  ~ 


|w"Ax|2 

w"Raw 


(40) 


Since  R  A'  is  only  an  approximation  of  the  inverse  of  RA,  the  denominator  of  (40) 
usually  cannot  be  further  simplified  as  has  been  done  in  (9). 


4.2  Two  Significant  Advantages  of  Proposed  ADPCA  over  STAP 

There  are  two  significant  advantages  of  the  proposed  ADPCA  over  STAP,  one  is 
computational  savings  and  the  other  is  its  robustness  when  reduced  samples  are  used 
in  parameter  estimation. 


4.2.1  Computational  Savings 

The  computational  savings  can  be  readily  estimated.  In  the  estimation,  we  only  count 
complex  multiplication/ division,  and  ignore  addition/ subtraction.  For  simplicity  we 
also  assume  that  no  special  digital  signal  processing  (DSP)  hardware  is  used,  so  that  all 
calculations  are  treated  with  the  same  weight.  Special  structures  of  matrices  such  as 
complex  conjugate  symmetry,  which  often  lead  to  special  algorithms  to  save 
operational  counts  (ops),  are  also  not  considered.  For  instance,  the  number  of  ops  for  a 
matrix  multiplication  Z  =  XY  is  Imn  where  Z  e  C"!X" ,  X  e  (C",x/  and  Y  e  C/x" . 

To  estimate  the  total  ops  required,  we  further  assume  that  the  fixed-window 
architecture  is  used  in  the  process.  The  fixed-window  architecture  only  uses  a  single 
covariance  matrix  estimated  from  a  section  of  K  range  samples  (secondary  data)  to 
process  a  section  of  K  range  cells  (primary  data).  In  contrast,  the  sliding- window 
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architecture  estimates  the  covariance  matrix  for  each  range  cell  (or  a  few  cells)  from  its 
surrounding  range  cells.  According  to  a  previous  study,  the  sliding-window 
processing  does  not  show  any  noticeable  SINR  improvement  compared  to  the  fixed- 
window  processing,  though  the  former  requires  a  much  more  considerable  amount  of 
computation  (Dong,  2006). 

Ops  required  for  STAP  and  the  proposed  ADPCA  algorithms  are  listed  in  Table  1  and 
Table  2,  respectively.  To  easy  understand,  calculations  are  divided  into  5  steps  (cl  to 
c5 )  for  STAP  and  9  steps  ( dx  to  d9 )  for  the  proposed  ADPCA.  Since  the  Doppler  of  the 
potential  target  signals  is  unknown,  a  search  in  Doppler  frequency  is  necessary.  In 
general  we  need  at  least  M  Doppler  filters.  The  total  ops  required  for  processing  K 
range  samples  using  the  fixed  window  architecture  for  STAP  is  therefore, 

c,  =  Cj  +  c2  +  M(c3  +  c4 )  +  KrgMcs  (41) 

Similarly,  the  total  ops  required  for  processing  K  range  samples  for  the  proposed 
ADPCA  is, 

dt  =dx+  d~,  +  d4  +  d  4  +  M  (r/5  +  )  +  Krg  (c/8  +  M  d 9)  (42) 


Table  1:  Ops  required  for  the  STAP  process. 


Algorithm 

Task 

Reference 

Eqn 

Ops 

Cj  -  form  R„  using  K  range 
samples 

1  K 

Ru=-Ix,xf  +SI 

K  k= 1 

(2) 

K(MN)2 

C2  -  compute  Ru ' 

r,;1 

(1) 

(MA)3 

C3  -  compute  Wopt 

wop,  =  R„e 

(1) 

(. MN )2 

C4  -  compute  A  ST4p0 

A- STAPO  =  Wopt^opt 

(9) 

MN 

C5  -  compute  Astap 

Iw"  xl2 

A  =  opt  1 

1  v  STAP  a 

^  STAPO 

(9) 

MN 
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Table  2:  Ops  required  for  the  proposed  ADPCA  process. 


Algorithm 

Task 

Reference 

Eqn 

Ops 

dl  -  form  A 

Vm-2  Tl[M~2  1 

L  m= 0  J  |_  m= 0  J 

(22) 

2M  aKN2a  +  2Na 

d 2  -  form  K  Ax's 

Ax  =  Bx 

(23) 

K(MaN2a) 

t/3  -form  Ra 

=^ZAx,Axf 

k= 1 

(27) 

K(MaNa)2 

d4  -  compute  Ra* 

Ra!  ~[diag(RAXr’ 

(28) 

MaNa 

d5  -form  U 

u  =  Be 

(25) 

MaN\ 

d6  -  compute  W 

w  =Rau 

(29) 

MaNa 

^7  '  ^ ADPCA0 

^-ADPCAO  =  W  ®-AW 

(30) 

C MANA)2+MANA 

ds  -  form  Ax 

Ax  =  Bx 

(23) 

MaN\ 

^9  '  ^  ADPCA 

wH Ax  2 

^  ADPCA  ~  * 

^  ADPCA  0 

(30) 

MaNa 

Note:  where  M  4  =  M  - 1  and  NA=N-\. 


Figure  3  shows  the  ratio  of  dt  /  ct  in  percentage,  as  a  function  of  the  number  of  pulses, 
M  ,  with  three  given  parameters,  the  number  of  antenna  elements  ( N  =  20,30,50 ),  the 
number  of  range  samples  ( K  =  100)  used  for  estimating  the  covariance  matrix  and  the 
number  of  range  cells  ( Kra  =500)  to  be  processed.  The  dominant  part  of  the  STAP 

algorithm  lies  in  the  computation  of  the  inverse  of  the  covariance  matrix,  which 
requires  ops  equal  to  the  cube  of  the  dimension  of  the  matrix.  On  the  other  hand,  the 
proposed  ADPCA  does  not  require  the  inversion  of  the  full  dimensional  covariance 
matrix,  so  its  ops  are  considerably  less.  In  general,  the  proposed  ADPCA  algorithm 
only  needs  about  5-10%  of  the  calculation  required  by  the  STAP  algorithm.  The  larger 
the  values  of  M  and  N ,  the  smaller  is  the  ratio  of  dt  /  c, . 

Furthermore,  as  mentioned,  the  inversion  of  the  covariance  matrix  dominates  the 
calculation  in  STAP.  The  method  of  the  complete  Gauss-Jordan  elimination  is  often 
used  to  numerically  compute  the  inverse  of  a  matrix  (Isaacson  and  Keller,  1966).  Since 
this  process  is  not  linear,  it  is  difficult  to  construct  dedicated  DPS  hardware  or  perform 
parallel  computation  to  implement  the  Gauss-Jordan  elimination.  On  the  other  hand,  it 
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can  be  seen  that  the  dominant  processes  for  ADPCA  are  linear.  A  linear  transform  can 
be  easily  implemented  using  either  general  DPS  hardware  or  parallel  computation.  In 
this  sense,  the  proposed  ADPCA  processing  should  have  an  even  bigger  computational 
advantage  than  the  simple  ops  estimation  and  comparison  shown  in  Figure  3. 


Figure  3:  Ratio  in  percentage  of  the  proposed  ADPCA  ops  to  the  STAP  ops  with  the  parameters 
of  K  =  100  and  Krg  =500. 

The  ops  required  for  the  Blum  ADPCA  processor  has  not  been  calculated.  Qualitatively 
its  computational  requirement  is  higher  than  that  of  the  proposed  ADPCA  but  much 
less  than  that  of  STAP.  It  will  be  seen  later  that  the  processor  does  not  perform  well,  so 
its  detection  loss  cannot  be  justified  by  its  savings  in  computation  requirements. 

4.2.2  Robustness  of  Proposed  ADPCA 

The  performance  of  any  adaptive  processor  depends  on  an  accurate  estimation  of  its 
parameters.  In  STAP  in  order  to  accurately  estimate  the  covariance  matrix,  twice  as 
many  iid  range  samples  as  the  size  of  the  covariance  matrix  itself  are  required  (Reed  et 
al,  1974,  Ward,  1994),  or  more  correctly,  about  twice  the  number  of  non-zero 
eigenvalues  of  the  matrix1.  Since  the  size  of  the  covariance  matrix  is  MN ,  a  significant 
number  of  range  samples  is  required  in  order  to  guarantee  a  satisfactory  performance 
of  STAP.  On  the  other  hand,  the  proposed  ADPCA  uses  a  covariance  matrix  whose 
size  is  only  N  - 1 .  This  suggests  that  the  number  of  range  samples  required  for  an 
accurate  estimate  of  the  parameter  matrix  should  be  significantly  smaller  than  that  for 
the  estimate  of  the  covariance  matrix  in  STAP.  A  physical  interpretation  for  the 
different  requirement  of  the  range  samples  is  that  STAP  only  utilises  averaging 
processing  in  the  fast-time  domain  (range  samples)  to  estimate  its  parameters  whereas 
the  proposed  ADPCA  employs  not  only  the  fast-time  but  also  the  slow-time  (pulse 


1Due  to  the  existence  of  thermal  noise,  the  non-zero  eigenvalue  means  that  the  eigenvalue  is 
greater  than  the  mean  power  of  the  thermal  noise. 
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samples)  averaging  processing  to  estimate  its  parameters.  Therefore,  in  the  case  of 
reduced  samples,  STAP  may  not  perform  due  to  a  poor  estimate  of  the  covariance 
matrix  while  the  ADPCA  still  robustly  performs.  More  detailed  analysis  is  given 
through  examples  in  the  next  section. 

4.2.3  Inherent  SINR  Loss  of  ADPCA 

We  need  to  point  out  that  the  proposed  ADPCA  is  a  dimensionality-loss  process.  Both 
the  number  of  effective  antenna  elements  and  the  number  of  effective  pulses  in  a  CPI 
are  reduced  by  one.  Therefore  theoretically  the  maximum  coherent  processing  gain  it 
can  achieve  is  101og10[(A-l)(M -1)]  dB  compared  to  the  gain  of  101og10(AM)  dB  for 
STAP. 


5.  Performance  Assessment 

A  dataset  generated  using  the  high  fidelity  airborne  radar  simulation  software, 
RLSTAP,  and  a  dataset  collected  by  the  MCARM  system  were  used  to  test  the 
performance  of  the  proposed  ADPCA.  The  results  of  the  Blum  ADPCA  were  also 
computed,  so  the  significant  improvement  of  the  proposed  processor  can  be 
appreciated.  The  results  of  STAP  serve  as  test  benchmarks  as  they  are  optimum. 

5.1  RLSTAP  Dataset 

Table  3  and  Table  4  list  parameters  of  the  radar,  platform  and  targets  used  in 
generating  the  simulated  dataset.  RLSTAP  calculates  clutter  returns  based  on  the 
United  States  Geological  Survey  Land  Use  Land  Cover  (USGS  LULC)  data  and  the 
Digital  Terrain  Elevation  (DTE)  data.  Figure  4  shows  the  LULC  data  superposed  onto 
the  DTE  data  for  the  Washington  D.C.  area  with  the  radar  beam  pattern.  Details  of  the 
clutter  model  used  in  RLSTAP  are  unknown. 
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Table  3:  Parameters  used  in  RLSTAP  for  generating  the  simulated  airborne  radar  dataset. 


Parameter 

Specification 

Radar 

Phased  array 

Linear  20-by-4  elements,  element  pattern  cos0'6  ( t/>) , 
azimuth  spacing  0.12m,  elevation  spacing  0.15m 
uniform  tapering  for  transmit 

Carrier  frequency 

1.2  GHz 

Polarisation 

W 

LFM  bandwidth 

2  MHz 

PRF 

2  KHz 

Number  of  pulses  per  CPI 

32 

Peak  power 

30kW 

Duty 

10% 

Sample  interval 

0.2  jus 

Look  direction 

30°  from  broadside  towards  nose,  horizontal 

Platform 

Height 

7  km 

Speed 

175  m/s 

Undesired  Signals 

Thermal  noise 

Gaussian 

Clutter 

Washington  D.  C.  area:  (38.1°N,  -76.9°E)  to  (39.6°N,  - 
75.2°E) 

Jamming 

None 

Table  4:  Parameters  of  targets. 


Parameters 

Target  1 

Target  2 

Target  3 

Height  (km) 

3 

3 

3 

Position  off  mainlobe  direction 

1  km  north 

0.5  km  south 

0 

Radial  velocity  (m/s) 

-100 

-20 

150 

Doppler  frequency  (Hz) 

-500 

860* 

-500 

Range  (km) 

50 

60 

70 

Range  bin  Number 

1948 

2256 

2598 

RCS  (sqm) 

1 

1 

1 

*Since  the  radar  looks  at  30°  from  broadside  towards  the  nose,  the  main  lobe  clutter  has  a 
Doppler  frequency  of  700  Hz.  Therefore,  target  2  is  the  most  difficult  to  detect  as  its  Doppler 
frequency  is  close  to  that  of  the  mainlobe  clutter. 


21 


DSTO-RR-0334 


Figure  4:  The  LULC  data  superposed  on  the  DTE  data  of  the  Washington  D.  C.  area. 

Before  showing  the  results,  it  is  worth  noting  that 

■  Due  to  the  effect  of  the  clutter  foldover,  clutter  returns  from  the  first  pulse 
statistically  differ  from  those  collected  from  the  following  pulses.  The  data  collected 
by  the  first  pulse  was  therefore  discarded  in  the  calculation.  The  number  of  pulses 
in  a  CPI  and  the  number  of  antenna  elements  in  the  azimuth  actually  used  in  the 
calculation  are  M  =  3 1  and  N  =  20,  respectively. 

■  According  to  the  parameters  given  in  Table  3,  (5  =  1 .46  so  the  simulation  was  under 
a  non-DPCA  condition. 

■  The  range  resolution  of  the  system  is  A r  =c/(2B)  where  c  is  the  speed  of  light  and 
B  the  bandwidth  of  the  LFM.  The  corresponding  range  sampling  interval  for  iid 
samples  should  be  At  =  V B  =  0.5/49 .  The  actual  range  sampling  interval  used  in 
RLSTAP  was  0.2/js  .  Therefore,  the  consecutive  range  samples  will  not  be  iid. 
Theoretically  every  third,  or  greater,  samples  should  be  iid  samples. 

■  The  DL-SMI  method  (Carlson,  1988)  was  used  in  the  STAP  analysis.  The  diagonally 
loaded  coefficient  was  set  to  be  50dB  below  the  mean  value  of  the  diagonal  elements 
of  the  covariance  matrix.  This  diagonal  loading  means  that  the  single-pulse  single¬ 
element  CNR  is  no  better  than  50dB  (the  single-pulse  single-element  CNR  for  the 
simulation  was  unknown  but  could  be  computed  from  the  data). 

■  Hamming  windows  were  used  for  both  the  spatial  and  temporal  steering  vectors 
e(/„)  and  e(ftd) in  the  processing. 
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Figure  5  and  Figure  6  compare  the  results  of  the  STAP,  DPCA  (non-adaptability), 
Blum's  ADPCA  and  the  proposed  ADPCA  obtained  from  using  267  non-consecutive 
range  samples  with  an  interval  of  3  (801:3:1600)2.  Figure  5  shows  the  signal  level  as  a 
function  of  both  the  range  bin  and  Doppler  bin,  whereas  Figure  6  is  a  view  of  Figure  5 
from  the  direction  perpendicular  to  the  range  bin,  i.e.,  the  Doppler  bins  are  collapsed 
onto  the  range.  It  was  found  that  STAP  produced  the  best  result  if  sufficient  range 
samples  were  used  in  forming  the  covariance  matrix.  It  can  be  seen  that  there  is  a  few 
dB  signal  loss  for  target  2  (the  mid  range  one,  that  is  the  most  difficult  to  be  detected, 
see  the  note  on  Table  4)  in  the  results  of  both  ADPCA  processors  (the  loss  of  the  Blum 
ADPCA  is  higher  though)  compared  to  the  result  of  STAP,  which  can  be  clearly 
viewed  from  Figure  6.  This  is  not  a  surprise  result,  as  STAP  is  optimal  if  the  covariance 
matrix  is  accurately  estimated.  The  DPCA  with  non-adaptability  processor  performs 
poorly  as  expected  since  the  radar  was  operated  under  a  non-DPCA  condition.  Further 
numerical  examples  will  not  include  this  DPCA  processor. 
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Figure  5:  Performance  comparison  between  (a)  STAP,  (b)  DPCA,  (c)  Blum's  ADPCA  and  (d) 
proposed  ADPCA  for  the  RLSTAP  data  using  267  range  samples  informing  the  covariance 
matrix. 


2The  use  of  800  (801:1600)  consecutive  range  samples  gave  a  very  similar  result. 
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Figure  6:  Performance  comparison  between  (a)  STAP,  (b)  DPCA,  (c)  Blum's  ADPCA  and  (d) 
proposed  ADPCA  for  the  RLSTAP  data  using267  range  samples  informing  the  covariance 
matrix.  The  results  are  plotted  in  the  signal  level  versus  range  with  the  Doppler  bins  collapsed 
onto  the  range. 


Reducing  the  number  of  range  samples  for  estimation  of  the  covariance  matrix  leads  to 
quite  different  results.  Figure  7  and  Figure  8  compare  the  results  of  STAP  and  the  two 
ADPCA  algorithms  using  100  non-consecutive  range  samples  (1301:3:1600).  It  can  be 
seen  that  the  clutter  level  has  lifted  significantly  for  the  STAP  due  to  insufficient  range 
samples  used  for  estimating  the  covariance  matrix.  Target  2  becomes  hardly  detectable 
for  the  Blum  ADPCA.  The  performance  of  the  proposed  ADPCA,  however,  remains 
almost  unchanged. 
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Figure  7:  Performance  comparison  between  (a)  STAP,  (b)  Blum's  ADPCA  and  (c)  proposed 
ADPCA  for  the  RLSTAP  data  using  100  range  samples  informing  the  covariance  matrix. 
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(a)  STAP  (b)  Blum's  ADPCA 
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Figure  8:  Performance  comparison  between  (a)  STAP,  (b)  Blum's  ADPCA  and  (c)  proposed 
ADPCA  for  the  RLSTAP  data  using  100  range  samples  informing  the  covariance  matrix.  The 
results  are  plotted  in  signal  level  versus  range  with  the  Doppler  bins  collapsed  onto  the  range. 

An  example  with  a  further  reduction  in  range  samples  is  shown  in  Figure  9  and  Figure 
10  in  which  only  50  (1451:3:1600)  non-consecutive  range  samples  were  used  in 
estimating  the  covariance  matrix.  It  can  be  seen  that  for  such  a  small  number  of  range 
samples  both  STAP  and  Blum's  ADPCA  stop  performing,  whereas  the  proposed 
ADPCA  still  robustly  works  with  little  degradation. 
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Figure  9:  Performance  comparison  between  (a)  STAP,  (b)  Blum's  ADPCA  and  (c)  proposed 
ADPCA  for  the  RLSTAP  data  using  50  range  samples  informing  the  covariance  matrix. 
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Figure  10:  Performance  comparison  between  (a)  STAP,  (b)  Blum's  ADPCA  and  (c)  proposed 
ADPCA  for  the  RLSTAP  data  using  50  range  samples  informing  the  covariance  matrix.  The 
residts  are  plotted  in  the  signal  level  versus  range  with  the  Doppler  bins  collapsed  onto  the 
range. 

In  Section  4  we  discussed  the  advantages  and  disadvantages  of  the  averaging 
processing  in  the  slow-time  domain.  In  order  to  see  the  effect  of  this  process  especially 
with  a  reduced  number  of  range  samples,  we  show  in  Figure  11  the  results  of  the 
proposed  ADPCA  without  slow-time  averaging,  i.e..  Mm) ,  m  =  (),•••,  M  -  2  are 
calculated  individually  using  (31).  Range  samples  used  are  the  same  as  those  of  Figure 
5,  Figure  7  and  Figure  9,  respectively.  It  can  be  seen  that  the  averaging  processing 
helps  significantly  for  the  case  of  reduced  range  samples. 
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Figure  11:  Results  of  the  proposed  ADPCA  without  using  slow-time  averaging  processing. 


It  is  known  that  in  order  to  accurately  estimate  the  covariance  matrix,  the  number  of 
samples  required  is  twice  the  dimension  of  the  covariance  matrix  itself  (Reed  et  al, 
1974,  Ward,  1994),  or  more  correctly,  about  twice  the  number  of  non-zero  eigenvalues 
of  the  matrix  (see  footnote  1).  Figure  12  shows  the  eigenvalues  of  the  covariance  matrix 
in  descending  order  (the  matrix  itself  was  formed  using  800  consecutive  range  samples 
(800:1600)  with  a  diagonal  loading  of  50dB  below  the  mean  value  of  the  diagonal 
elements  of  the  covariance  matrix).  The  number  of  eigenvalues  above  the  noise  level  is 
about  200.  Therefore,  theoretically  STAP  requires  about  400  iid  range  samples  in  order 
to  accurately  estimate  the  covariance  matrix,  or  about  250  iid  range  samples  if  we  treat 
those  eigenvalues  that  are  60dB  below  the  maximum  eigenvalue  to  be  noise. 
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Figure  12:  Eigenvalues,  in  descending  order,  of  the  sampled  clutter  covariance  matrix  using 
800  range  samples. 

The  proposed  ADPCA,  on  the  other  hand,  uses  matrices  whose  dimension  is  N  —  1  to 
estimate  its  parameters.  As  a  result,  it  requires  significantly  fewer  iid  range  samples.  In 
addition,  STAP  only  uses  the  fast-time  (range)  averaging  processing  to  estimate  its 
parameters,  whereas  the  proposed  ADPCA  employs  not  only  the  fast-time  but  also  the 
slow-time  (pulse)  averaging  processing  to  estimate  its  parameters.  Therefore,  the 
proposed  ADPCA  is  more  robust  in  the  case  of  reduced  samples. 

5.1.1  Comparison  with  PAMF 

The  proposed  ADPCA  algorithm  looks  similar  to  the  PAMF  algorithm  (Roman,  et  al, 
2000,  Michels  et  al,  2003,  Dong,  2006).  However,  there  exist  differences.  First,  the 
proposed  ADPCA  seems  to  be  equivalent  to  PAMF  with  filtering  order  equal  to  one. 
PAMF  with  the  filtering  order  equal  to  one  usually  does  not  perform  very  well  but 
ADPCA  does.  The  key  point  is  that  ADPCA  utilises  the  principle  of  DPCA,  whereas 
PAMF  usually  has  to  use  a  higher  order  of  filtering  to  gain  a  good  performance. 
Typical  filtering  orders  are  3  to  5  (Roman,  et  al,  2000,  Dong,  2006).  An  increase  in  the 
filtering  order,  however,  also  increases  the  size  of  the  matrix  to  be  inverted  (the  size  of 
the  matrix  to  be  inverted  is  the  filter  order  times  the  number  of  antenna  elements), 
which  in  turn  increases  both  the  computational  cost  and  the  size  of  sample  data.  PAMF 
requires  the  inversion  of  diagonal  block  matrices  of  size  N  whereas  the  ADPCA  does 
not  have  such  requirement.  Figure  13  and  Figure  14  show  the  performance  of  PAMF 
with  different  filtering  orders  in  which  100  (1301:3:1600)  non-consecutive  range 
samples  were  used.  It  can  be  seen  that  when  the  filtering  order  reduces  to  one,  which  is 
equivalent  to  the  ADPCA  algorithm,  the  processor  fails  to  detect  target  2.  However, 
there  is  little  SINR  loss  for  target  2  if  the  filter  order  increases  to  5. 
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Figure  13:  Performance  of  the  PAMF  processor  with  different  filtering  orders. 
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(a)  Filtering  order  =  5 
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(b)  Filtering  order  =  1 


Figure  14:  Performance  of  the  PAMF  processor  is  dependent  on  the  filtering  order.  The  results 
are  plotted  in  the  signal  level  versus  range  with  the  Doppler  bins  collapsed  onto  the  range. 


5.2  MCARM  Dataset 

Detailed  descriptions  of  the  MCARM  system  can  be  found  elsewhere  (Sloper  et  al, 
1996,  Fenner  and  Hoover,  1996).  Some  of  the  MCARM  data  analyses  are  also  available 
(MITRE,  1999,  RAFDCI,  1999,  Sarker  et  al,  2001).  The  dataset  #5-575  collected  by  the 
MCARM  system  was  used  in  this  report.  The  radar  and  platform  parameters  of  #5-575 
are  given  in  Table  5  and  Table  6,  respectively. 


Table  5:  MCARM  radar  parameters. 


Frequency/ 

Polarisation 

PRF 

CPI 

Pulse 

Width 

Duty 

Cycle 

Range 

Resolution 

PRI 

(irs  /gates) 

1240  MHz/VV 

1984 

128 

50.4  |rs 

10% 

0.8  |rs 

504/630 

Table  6:  MCARM  platform  parameters. 


Height 

Velocity 

Illumination 

Crab  Angle 

Antenna  Tilt  Angle 
from  Horizontal^ 

3488  m 

100.1  m/s 

Side-looking 

7.28° 

5.11° 

The  receiver  of  the  MCARM  system  consists  of  22  receiving  channels  (modules) 
organised  in  two  azimuthally  identical  rows.  We  treated  data  received  by  one  row  as 


3  This  angle  is  the  sum  of  the  antenna  tilt  angle  relative  to  the  platform  plus  the  recorded 
platform  roll  angle. 
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secondary  data  and  the  data  received  by  the  other  as  primary  data.  The  physical 
position  difference  between  the  two  rows  will  in  general  induce  a  phase  difference 
between  signal  arrivals  at  the  two  rows.  This  however  should  impose  no  effect  as  far  as 
generating  the  covariance  matrix  is  concerned.  We  also  only  used  the  data  collected  by 
the  first  50  pulses  in  the  process,  hence,  =  11  and  M  =  50  in  this  MCARM  data 
study. 

The  antenna  element  horizontal  spacing  is  0.1092  m.  With  the  parameters  given  in 
Table  5  and  Table  6,  we  found  J3  =  0.924 ,  indicating  that  the  data  was  collected  under 
a  non-DPCA  condition. 

Figure  15  shows  the  clutter  profile  in  range  after  the  1  /  R}  range  effect  is  compensated. 
The  area  illuminated  by  the  radar  mainlobe  is  mainly  farmland  with  scattered  houses 
(range  bins  200-400  and  500-630,  for  example)  and  bay  water  (range  bins  400-500,  for 
example).  The  signal  in  range  bin  68  is  a  replica  of  the  transmitter  signal  and  the  useful 
clutter  data  are  from  around  range  bin  200  and  beyond. 


Figure  15:  Clutter  profile  in  range  after  the  1/R3  range  effect  was  compensated. 

A  moving  target  in  range  bin  299  (a  range  ambiguous  target)  has  been  detected  and 
reported  previously  (Dong,  2005).  Range  cells  289  to  309  were,  therefore,  excluded  in 
the  secondary  data.  An  artificial  moving  target  signal  was  also  injected  in  the  primary 
data  in  range  bin  500  with  a  constant  amplitude  30dB  below  the  mean  clutter  of  that 
bin  and  a  Doppler  frequency  of  -200  Hz. 


To  evaluate  the  performance  of  the  proposed  ADPCA  processor,  we  first  formed  the 
covariance  matrix  using  nearly  all  the  possible  secondary  data,  range  bins  200  to  600 
with  the  exclusion  of  bins  289  to  309.  The  results  are  shown  in  Figure  16.  It  can  be  seen 
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that  all  three  processors  have  successfully  detected  both  the  genuine  target  in  range  bin 
299  and  the  artificial  target  in  range  bin  500.  However,  the  detection  of  the  artificial 
target  in  the  STAP  processing  as  well  as  the  Blum  ADPCA  suffers  a  few  dB  signal  loss 
compared  to  the  detection  of  the  ADPCA  processing.  This  can  be  more  clearly  viewed 
in  Figure  17  where  the  figures  are  plotted  in  the  signal  level  versus  range  with  the 
Doppler  bins  collapsed  into  the  range.  Other  evidence  readily  seen  from  Figure  16  is 
that  there  are  higher  signal  levels  in  the  mainlobe  clutter  Doppler  vicinity  (note  that 
due  to  the  aircraft's  crab  angle,  the  Doppler  of  the  mainlobe  clutter  is  about  150  Hz)  in 
the  results  of  both  the  STAP  and  Blum's  ADPCA,  indicating  that  the  clutter  is  not 
sufficiently  suppressed,  possibly  due  to  the  insufficient  range  samples  used  for 
covariance  matrix  estimation,  although  all  possible  range  samples  have  been  used.  No 
such  evidence  however  is  seen  in  the  result  of  the  proposed  ADPCA. 
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Figure  16:  Detection  results  ofMCARM  data  by  the  use  of  (a)  STAP,  (b)  Blum's  ADPCA  and 
(c)  proposed  ADPCA.  The  covariance  matrix  was  formed  using  approximately  400  range 
samples. 


34 


DSTO-RR-0334 


Range  bin  RanSe  bm 

(a)  STAP  (b)  Blum's  ADPCA 


Range  bin 

(c)  Proposed  ADPCA 


Figure  17:  Detection  results  ofMCARM  data  by  the  use  of  (a)  STAP,  (b)  Blum's  ADPCA  and 
(c)  proposed  ADPCA.  The  covariance  matrix  was  formed  using  approximately  400  range 
samples.  The  results  are  plotted  in  the  signal  level  versus  range  with  the  Doppler  bins  collapsed 
into  the  range. 

To  test  the  advantages  of  the  proposed  ADPCA  in  the  case  of  reduced  samples,  we 
greatly  reduced  the  number  of  samples.  The  LFM  bandwidth  of  the  MCARM  system 
was  1  MHz  and  the  actual  range  sample  interval  was  0.8  ps.  To  ensure  that  samples 
were  statistically  independent,  we  used  every  second  range  sample  as  the  iid  samples. 
Results  with  the  use  of  21  (200:2:240)  non-consecutive  range  samples  to  form  the 
covariance  matrix  are  shown  in  Figure  18.  Comparing  Figure  18  with  Figure  16,  we  can 
see  that  the  output  of  STAP  has  been  seriously  distorted  due  to  reduced  range  samples 
which  cannot  provide  an  accurate  estimate  of  the  covariance  matrix  and  hence 
suppress  the  mainlobe  clutter.  The  situation  of  the  Blum  ADPCA  is  similar,  but  the 
distortion  is  less  severe.  On  the  other  hand,  the  differences  between  results  of  the 
proposed  ADPCA  shown  in  Figure  16  and  Figure  18  are  minor,  indicating  that  the 
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parameters  of  ADPCA  can  be  estimated  at  a  reasonably  accurate  level  even  using 
reduced  samples.  The  plots  of  the  signal  level  versus  range  with  the  Doppler  bins 
collapsed  into  the  range  are  shown  in  Figure  19. 
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Figure  18:  Detection  results  ofMCARM  data  by  the  use  of  (a)  STAP ,  (b)  Blum's  ADPCA  and 
(c)  proposed  ADPCA.  The  covariance  matrix  was  formed  using  21  range  samples. 
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Figure  19:  Detection  results  ofMCARM  data  by  the  use  of  (a)  STAP,  (b)  Blum's  ADPCA  and 
(c)  proposed  ADPCA.  The  covariance  matrix  was  formed  using  21  range  samples.  The  residts 
are  plotted  in  the  signal  level  versus  range  with  the  Doppler  bins  collapsed  into  the  range. 

Like  STAP,  the  proposed  ADPCA  is  fully  adaptive  to  changes  in  the  clutter 
environment  and  system  parameters  provided  the  effects  of  these  changes  have  been 
included  in  the  sample  data  used  for  the  parameter  estimation.  In  the  above  MCARM 
data  analysis,  we  have  seen  that  the  proposed  ADPCA  is  able  to  cope  with  all  effects 
such  as  the  antenna  pattern  distortion/  degradation  unavoidably  induced  by  the 
aircraft,  the  shift  of  the  Doppler  frequency  of  the  mainlobe  clutter  induced  by  the 
crabbing  angle  (as  large  as  7°  for  the  MCARM  dataset),  the  clutter  notch  widening 
induced  by  the  clutter  intrinsic  motion  and  range  foldover  (see  the  report  of  Dong, 
2005,  for  detailed  discussion  of  the  effect  of  the  range  foldover)  as  well  as  various  other 
spatial  and  temporal  decorrelation  in  the  MCARM  system.  However  we  should  not 
confuse  the  ability  to  cope  with  changes  and  the  degradation  of  the  radar  performance 
induced  by  these  changes.  The  latter  is  the  inherent  result  of  the  effects  and  cannot  be 
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recovered  by  any  means.  For  instance,  a  crab  angle  of  7°  of  the  MCARM  system  causes 
the  shift  of  the  mainlobe  clutter  Doppler  to  be  about  150  Hz.  As  a  result,  any  moving 
targets  whose  Doppler  frequencies  are  in  the  vicinity  of  150  Hz  and  whose  echoes  are 
not  as  strong  as  the  mainlobe  clutter  will  become  undetectable.  The  clutter  notch 
widening  caused  by  the  various  decorrelation  effects,  clutter  intrinsic  motion  and 
range  foldover  will  also  reduce  the  radar's  effective  detection  region. 


6.  Conclusions 

In  this  report  we  have  proposed  an  algorithm,  the  adaptive  displaced  phase  centre 
antenna  (ADPCA),  to  process  airborne  phased  array  radar  data  for  moving  target 
detection.  Its  parameters  are  adaptively  calculated  by  the  autoregressive  (AR)  process. 
Although  the  name  is  not  new,  this  proposed  ADPCA  is  fundamentally  different  from 
the  existing  ADPCA  as  the  former  is  an  optimum  processor  whereas  the  latter  is  not. 
The  proposed  algorithm  not  only  significantly  reduces  the  number  of  samples  required 
for  estimating  its  parameters  but  also  demands  dramatically  less  computational  effort 
with  little  sacrifice  of  the  SINR  loss  in  comparison  to  the  conventional  space-time 
adaptive  processing  (STAP).  The  algorithm  does  not  have  any  assumptions,  so  it  is 
fully-adaptive.  The  performance  of  the  algorithm  has  been  assessed  using  two  airborne 
radar  datasets,  one  generated  using  the  high  fidelity  airborne  radar  simulation 
software,  RLSTAP,  and  the  other  collected  by  the  MCARM  system.  These  two  datasets 
were  carefully  chosen  to  cover  various  issues.  First  of  all,  the  datasets  did  not  satisfy 
the  DPCA  condition,  so  that  the  adaptability  of  the  proposed  ADPCA  could  be 
examined.  Secondly  the  radar  did  not  look  in  the  broadside  direction  in  the  RLSTAP 
dataset.  The  second  dataset  was  collected  from  a  flight  trial  of  MCARM  and  included 
effects  of  aircraft  crabbing  motion.  Other  decorrelation  effects  including  radar 
instability,  clutter  intrinsic  motion,  range  foldover  and  interference  caused  by  the 
aircraft  etc.  were  hence  automatically  included  in  the  dataset.  Therefore  the  evaluation 
of  the  proposed  processor  is  realistic.  The  results  of  STAP  have  served  as  benchmarks 
in  the  evaluation. 

It  has  been  found  that  the  proposed  ADPCA  performs  nearly  as  well  as  STAP, 
suffering  at  most  a  few  dB  of  processing  gain  loss  in  the  vicinity  of  the  Doppler  of  the 
mainlobe  clutter.  However  if  there  are  insufficient  clutter  samples  to  accurately 
estimate  the  covariance  matrix,  the  performance  of  STAP  is  severely  degraded  whereas 
the  proposed  ADPCA  still  performs  as  well  as  before.  Mathematically  STAP  requires 
estimation  of  a  covariance  matrix  whose  size  is  the  product  of  the  number  of  antenna 
elements  and  the  number  of  pulses  in  a  CPI.  On  the  other  hand,  the  proposed  ADPCA 
estimates  its  parameters  by  the  use  of  a  covariance  matrix  whose  size  is  no  more  than 
the  number  of  antenna  elements.  In  addition,  in  estimating  parameters,  STAP  uses 
averaging  processing  only  in  the  fast-time  domain,  whereas  ADPCA  utilises  averaging 
processing  in  both  the  fast-time  and  the  slow-time  domains.  These  two  differences 
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make  the  proposed  ADPCA  more  robust  and  require  far  less  sample  data.  In  general, 
the  parameters  of  the  ADPCA  can  be  satisfactorily  estimated  once  the  number  of  range 
samples  is  equal  to  or  greater  than  twice  the  antenna  elements.  The  Blum  ADPCA,  on 
the  other  hand,  is  not  an  optimum  processor  and  its  performance  usually  is  poor  and 
suffers  significantly  especially  when  the  target's  Doppler  is  close  to  that  of  the 
mainlobe  clutter. 

The  operational  counts  (ops)  required  for  the  proposed  ADPCA  has  been  estimated.  In 
general,  it  only  requires  5-10%  of  the  computation  of  STAP.  Since  most  of  the 
computation  for  the  ADPCA  algorithm  is  linear  transforms,  parallel  processing  and /  or 
hardware  realisation  can  be  easily  implemented.  In  contrast,  the  dominant  calculation 
of  the  STAP  algorithm  is  the  inversion  of  the  covariance  matrix  which  is  not  a  linear 
transform  and  limits  the  application  of  parallel  processing.  In  this  sense,  computational 
savings  of  the  ADPCA  algorithm  are  even  greater  than  the  simple  ops  estimation  and 
comparison  presented  in  the  report. 
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